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Prediction of a surface magnetic moment in a-uranium 
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Recently, there has been an increased interest in first-principles calculations of the actinides as 
well as in finding the new materials which display surface magnetism. We predict the existence of a 
magnetic moment on the uranium (001) surface by performing density functional calculations for a 
slab geometry in the generalized gradient and local spin density approximations with included spin 
orbit coupling. The ferromagnetic phase is energetically favored for all geometries. The calculated 
total magnetic moment, 0.65/^s, is stable on films of different thickness and it should be observable 
experimentally. 



PACS numbers: 75.70.Ak, 75.70.Rf, 71.20.-b, 71. 15. Mb 



The beginnings of the field of surface magnetism can 
be traced back to the end of the 1960s with the claims of 
magnetic dead layers of some transition elements. In sub- 
sequent years it has been shown both experimentally and 
by first-principles calculations that Fe(OOl), Ni(OOl) and 
Cr(OOl) have considerable enhancement of magnetic mo- 
ments at the surface with respect to the bulk value, 0. 
As a rule one expects the tendency toward magnetism 
to be increased near metal surfaces, because of the nar- 
rowing of the density of states which yields a Stoner en- 
hancement in the susceptibility. This lead to the search 
for surface magnetism in materials which are nonmag- 
netic in the bulk. There have been many attempts to 
find such materials (vanadium, for example), but most 
of them proved unsuccessful. To our knowledge the only 
experimentally clear case of a magnetic surface on a non- 
magnetic metal is Rh(lOO), In contrast, there are 
many examples of predicted transition metal monolayers 
on noble-metal substrates such as: Ti, V, Tc, Ru, and 
Rh on Ag or Au, [3|. 

Uranium would appear to be a possible candidate for 
surface magnetism since there is a rather large change 
in density of states from surface to bulk, as pointed out 
by Hao et al. 0. Also, a comparatively modest expan- 
sion of the lattice leads to a spontaneous moment in 
the bulk which we have verified, [^. In addition, it is 
near the transition (which occurs beyond plutonium) in 
which the 5f electrons become localized, 0- We also note 
that many U alloys are magnetic, indicating the possible 
presence of a nearby magnetic instability. For example, 
URhAl m and UPtAl j3l are ferromagnetic with spin 
moments of order l/is. In these structures the U atoms 
are in hexagonal planes (which include the transition el- 
ement) and are separated by planes with no U. 

We have explored the possibility for magnetism in ura- 
nium within the framework of the density functional the- 
ory (DFT). For uranium and the lighter actinides, DFT 
has been successful at predicting structural properties 



[3111113 ■ For plutonium (and presumably for the heav- 
ier actinides) calculations which go beyond DFT (such as 
dynamical mean field theory) are required 13]. The cal- 
culations have been done using the orthorhombic struc- 
ture (a-uranium) , Fig^ This is the stable structure be- 
tween 43K and 940K. At higher temperatures, it passes 
through tetragonal and then bcc structures before melt- 
ing at 1408K. Below 43K the lattice undergoes a charge 
density wave (CDW) distortion [3. The CDW distor- 
tions are only «0.027 A, Q, and were neglected. The 
orthorhombic structure may be viewed as a distorted bcc 
structure with the (001) plane similar to bcc (110) and 
with alternate (110) planes shifted along the b axis by 
(1/2 — 2y)b. A similar transformation was described by 
Axe et al. in The bcc structure would have y =1/4, 
a = Ofccc, and b = c = \/2ai,cc- At the same volume 
af)cc=3.405 A. The orthorhombic structure has a =2.836 
A, b =5.866 A, c =4.935 A, and y =0.1017, It be- 
longs to the non-symmorphic space group Cmcm. We 
confirmed previous calculations, 11], which show that 
this structure is lower in energy than bcc by 0.22 eV. 
We have also found that the bcc structure is mechani- 
cally unstable in the sense that the bcc energy is a rel- 
ative maximum with respect to tetragonal distortions (a 
pho non calculation would lead to negative frequencies), 
[l7|. For the calculations reported here we use the unre- 
laxed experimental geometry. 

In this pap er we solve the DFT equations using the 
WIEN2k, [l8j, implementation of the full potential lin- 
ear augmented plane wave (FLAPW) method, J^], in a 
supercell geometry. To simulate a surface, we varied the 
length of the vacuum layer until no changes were noticed 
in energy and density of states and we kept it at 18 Bohr. 
Our free-standing monolayer calculations were done us- 
ing 2 atoms per unit cell on a (001) surface, experimental 
lattice constants, 150 k points in the irreducible Brillouin 
zone and energy cutoff 12.5 Ry. For the film calculations 
(1 atom per surface cell), number of k-points varied from 
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FIG. 1: a -uranium structure. 



306 for the three-layer film to 153 for the seven-layer film. 
The muffin tin radius was kept fixed at 2.55 Bohr. For 
the exchange and correlation we used both the general- 
ized gradient (GGA) , '2^ , and the local spin density ap- 
proximations (LSDA) for the monolayer and only GGA 
for the film calculations. In the actinides it is important 
to include several corrections to the usual scalar relativis- 
tic formulation of the FLAPW method. The spin orbit 
splitting of the 5f level in the isolated atom is 0.82 eV. 
This is comparable, though smaller than the f bandwidth 
which is of order 2 eV. Hence it is important to include 
spin orbit effects explicitly. This was done treating the 
spin orbit coupling in a second variation, Our bulk 
energy bands and density of states are in good agree- 
ment with fully relativistic results, In addition, the 
semicore 6s and 6p levels have bandwidths in the solid 
of order 2 eV so it is important to include them in the 
basis. In the Wien code they are treated by adding local 
orbitals at the appropriate energies, [2l| . 





GGA 


GGA -f-SOC 


LSDA 


LSDA -\- SOC 


FM 


-0.333 


-0.226 


-0.241 


-0.163 


AFM 


-0.206 


-0.136 


-0.141 


-0.089 



TABLE L Calculated energy differences from the paramag- 
netic phase for the ferromagnetic (FM) and antiferromagnetic 
(AFM) state with and without spin-orbit coupling in GGA 
and LSDA for a monolayer in eV/atom. 



We find that the surface does indeed support a mag- 
netic moment. Our results for a monolayer, presented in 
Table U demonstrate that the ferromagnetic (FM) phase 
is favored over antiferromagnetic (AFM) and paramag- 
netic (PM) phases. For the GGA with spin-orbit coupling 
(SOC) included, it hes 0.23 eV below the PM phase and 
0.09 eV below the AFM phase. This result qualitatively 
does not depend on the exchange-correlation potential 



GGA GGA+SOC LSDA LSDA -f SOC 


FM - 


sphere 


2.42 2.11 (1.27) 


2.13 


1.97 (1.13) 


FM - 


interstitial 


0.85 0.74 


0.71 


0.62 


FM - 


total per atom 


3.27 2.01 


2.84 


1.75 


AFM 


- total per atom 


2.25 1.92 (0.79) 


1.89 


1.69 (0.59) 



TABLE II: Calculated spin (spin -|- orbital) magnetic mo- 
ments on spheres and in the interstitial region for the FM 
and AFM state with and without spin-orbit coupling in GGA 
and LSDA for a monolayer in units oi \ib- Interstitial contri- 
bution is given per atom. 



and inclusion of SOC. The direction of magnetization in 
the case with SOC is [001], perpendicular to the surface. 
The values of spin and total (orbital moment taken into 
account) magnetic moments on spheres and in the inter- 
stitial region for the FM and AFM phase are given in 
Table Ull (Interstitial contributions come out as a conse- 
quence of the LAPW method's partitioning of space into 
atomic spheres and the interstitial region). We can see 
that the total moments per atom on the uranium mono- 
layer are surprisingly large, 2.01 and 1.75 [Ib in GGA and 
LSDA, respectively. Total moment per atom is obtained 
by adding spin and orbital moments from a sphere to the 
spin interstitial moment. The orbital moments in the in- 
terstitial region vanish. We also compare the magnitude 
of the magnetic moments for the case without SOC. The 
moments are significantly larger in that case, similar to 
4(i and hd (Ru, Rt, Pd, Os, Ir, Pt. . .) mono and double 
layers on Ag(OOl) and Au(OOl), H^. 

Magnetic moments are similar in GGA and LSDA 
approximations. The influence of the type of the ex- 
change and correlation approximation on the calculations 
of magnetic effects on surfaces is a subject of a discussion. 
For example, for chromium (001), |23, it was found that 
GGA tends to overestimate magnetic moments in bulk 
and at surface. However, the use of the LSDA for the de- 
scription of a magnetic state also gave debatable results 
(bulk iron, for instance). All results presented further in 
this work have been obtained using the GGA. 

The results for thick films are given in Table UTTI Since 
the calculated surface moment is stable upon the change 
in the slab thickness, we have performed calculations for 
up to seven layers. The constant value of the magnetic 
moment for the three cases is 0.65/iB. We show in Ta- 
ble mil the spin and orbital moments in the outermost 
sphere and in the interstitial region for films of varying 
thickness. Interstitial moments are constant for the films 
of different thickness, from which we conclude that they 
are mostly at the surface. They are given per total calcu- 
lational cell, so dividing them by two (two surfaces) and 
adding to the sum of spin and orbital moment, we get 
the total magnetic moment per atom. The moments di- 
minish rapidly as a function of distance from the surface, 
so that, for example, for the five-layer film the total mo- 
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ments on the spheres are: 0.46, -0.03, and -0.04/1^. The 
orbital moments are relatively large. This is consistent 
with the findings for uranimn compomids, 0. In Table 
liVl we present the energy differences from the PM to the 
FM phase for different films. 





1 layer 


3 layers 


5 layers 


7 layers 


spin moment 


2.11 


0.84 


0.84 


0.84 


orbital moment 


-0.38 


-0.38 


-0.38 


-0.38 


interstitial moment 


0.74 


0.39 


0.38 


0.38 


total moment 


2.01 


0.66 


0.65 


0.65 



TABLE III: Spin and orbital moments on spheres and in the 
interstitial region at the outermost layer for slabs of different 
thickness, with spin-orbit coupling, given in hb- For 3-7 
layers the total moment is taken to be spin -I- orbital -1-1/2 
interstitial. 



26. SA"^ (w30% increase). The moments increase with the 
further expansion. 

In Fig|3| we show the spin-polarized density of states 
for a monolayer, which has a very sharp difference be- 
tween spin up and down states. Figure 01 shows the spin- 
polarized density of states, for the surface layer of the 
7-layer film. 

Experimentally, the surface magnetism can be detected 
in a variety of ways including, for example, spin-polarized 
photoemission [3| , linear or circular dichroism j3| and sur- 
face magneto-optic Kerr effect [Isj . Given the size of our 
computed moments, these effects should be readily ob- 
servable. 

N. S. would like to thank E. Vescovo for many inter- 
esting and helpful discussions. Supported by US De- 
partment of Energy under Contract No. DE-AC02- 
98CH10886. 



1 layer 3 layers 5 layers 7 layers 
-0.226 -0.011 -0.016 -0.020 



TABLE IV: Energy differences from the paramagnetic state 
for different slab thickness geometries in eV/atom. 

In order to justify our use of the unrelaxed experimen- 
tal geometry in this work, we contracted the surface layer 
of a 5-layer slab. The energy minimum occurs at the con- 
traction of 2.66%. We used the muffin tin radius of 2.35 
Bohr and 72 k points in the irreducible Brillouin zone. 
The spin moment on a sphere of the relaxed structure 
was reduced slightly to 0.65/.t_B. 

We also calculated the surface energy using the un- 
relaxed geometry and including the SOC. It was found 
to be 0.12 eV/A2. This is in a resonable agreement 
with the only available experimental result of 0.094eV/A^ 
|23 |. Previous calculations of KoUar et ai, [l^, gave 0.16 
eV/A2. 

The existence of the surface magnetic moment can be 
ascribed to the reduction of the coordination number 
which causes narrowing and enhancing of the density of 
states at the Fermi level. This increase in the local den- 
sity of states, nioc(Ep), results in fulfilling the Stoner 
criterion nioc{Ep)I > 1, (/is the exchange integral). For 
reasonable values of I (~ 0.2 - 0.4 eV), |22|, this change 
in nioc is sufficient to drive the system magnetic, as we 
observe. This can be seen from the density of states for a 
paramagnetic 7-layer film and bulk shown in Fig. [21 The 
local density of states at the Fermi level changes from 
2.45 for bulk to 4.01 states/eV for the surface layer. It is 
even more pronounced for a monolayer, where PM local 
density of states at Fermi level reaches the value of 11.38 
states/eV, Fig. O A similar effect of an emerging mag- 
netic moment can be produced by a uniform expansion 
of the original lattice. For example, we have found that 
there is a magnetic moment of O.l/is at the volume of 
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FIG. 2: Density of states for paramagnetic 7-layer film and 
bulk. 
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FIG. 3: Density of states for paramagnetic and magnetic 
monolayer. 
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FIG. 4: Spin-polarized density of states for the outermost 
layer of the 7-layer film. The solid blue line represents the 
spin down density of states and the dot-dashed red line spin 
up density od states. 
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